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PART TWO 


GRAPHICAL ECLIPSE DATA 
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PRECEDING PAGE BEANIK EILW 


This section contains the eclipse data in graphical form. On 
facing pages are a work sheet (left side) and the final eclipse light 
curve (right side). (Figures 2-1 and 2-2 illustrate.) These figures 
also contain a number of explanatory comments which pertain to all 
the figures in Part Two. Figure 2-1 illustrates the "curve-fit" 
technique of removing the sky background from the uncorrected object 
counts after opening the telescope to 200-inch aperture. Prior to 
opening, the telescope aperture to 200-inches, the sky background was 
removed from the uncorrected object counts by subtracting the weighted 
sky count (as observed in the sky aperture of the spectrometer) from 
the object count (The "dual aperture" technique). The sky count was 
v;eighted by the balance factor appropriate to the particular channel. 
The balance factor was determined from the post disappearance data 
in the method discussed in Part One. 

In the case of the February ingress and March egress data, the 
dual aperture technique was used to correct the data after opening 
to 200-inch aperture as well. The curves are ordered by date-- 
February 8, March 10, April 6 and then by wavelength from blue to 
red. 

Error bars are shown on the data from the time the telescope was 
first opened to 200-inch aperture until the la level above zero 
counts. The data points are connected from the beginning until the 
data first falls below la. A discussion of errors is found in Part 
One, Section VI. 
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Figure 2-2 Annotated Eclipse Light Curve 
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TIME ORIGIN, FEBRURRY 8, 1971 
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TIME ORIGIN, FEBRUARY 8. 1971 13 HR 39 MIN 
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time origin. FEBRUARY 8. 1971 
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time origin, march 10, 1971 10 HR 56 MIN tU.T 
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ABSTRACT 

Observations of five Jovian satellite eclipses were obtained 
during February, March, and April 1971 with the 200-inch Hale telescope 
and multichannel spectrometer. Eclipse light curves in 20 wavelength 
bandpasses were obtained for each eclipse. The higher counting rate 
channels were statisitically significant during 10 magnitudes of 
darkening. Well developed refractive tails were evident in the April 6 
eclipse of Europe. 

Preliminary analysis of these data has been performed. The 
eclipse radius of Europa was found to be 1200 - 100 km. Employing a 
Jovian scale height of 10 km, extinction due to an aerosol component 
homogeneously mixed with the gas above the cloud tops, and total extinc 
tion due to cloud tops at the 7.58 x 10^^ particle/cc level it was 
possible to theoretically duplicate the observed eclipse light curves 
in wavelength bandpasses not characterized by molecular band absorption 
A wavelength variation in the aerosol extinction factor was found. 

Observations in several channels also manifested the presence of 
molecular bands. Strong NH3 and CH^ absorption was noted in several 
channels. Evidence is presented supporting the tentative conclusion 
that the 83, fi rst-order-forbidden , pressure induced dipole line at 
A 6420 ^ has been observed. These data will enable new determinations 
of the abundances of several of the atmospheric constitutents. 
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COVER 

The cover photograph shows Gan^miede during an eclipse ingress. 
The trailing limb of the satellite is just crossing the 
pen umbra -umbra boundary. After this configuration, the 
satellite is illuminated only by refracted light. Ganymede 
may also be observed during eclipse egress prior to occultation. 
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PART ONE 


I. INTRODUCTION 

Telescopic observations of the planet Jupiter at visual and near 
infrared wavelengths reveal a series of banded structures. These 
belts and zones vary in brightness with zenographic longitude and 
latitude, with time and with wavelength. In addition, they exhibit 
fine structure. At times the surface clouds show a pattern or regu- 
larity in their variations and multiple velocities have been observed 
indicating a shear velocity within the visible surface. These data 
point to a diversity of composition and structure within the visible 
1 ayers . 

Spectrograph! c observations of Jupiter have revealed the existence 
of minor amounts of Methane and Ammonia. In addition, the quadrupole 
rotation-vibrational lines of molecular hydrogen confirm the existence 
of ^ 2 - Further species have not yet been detected spectroscopically 
even though there is strong reason to expect other gasses. The earth's 
atmosphere makes observations of these constituents difficult. Either 
they are hidden by telluric lines and bands or the features are too 
weak to be detected by normal incidence observations. 

Other methods for determining the composition and physical structure 
of the atmosphere of Jupiter include observations of spacecraft and 
stellar occultations and observation of eclipses of the Galilean satel- 
lites. Flyby orbiter and direct entry probes could, of course, also be 
employed. The Pioneer and Outer Planet Exploration programs of NASA 



are pursuing these approaches. The ff-Arietis and the recent /3-Scorpius 
occultations by Jupiter illustrate the occultation approach. 

The more frequent phenomena are the Galilean Satellite Eclipses. 
While traditional spectroscopic measurements of the Jovian atmosphere 
at normal incidence reach to the cloud tops and occultations sample the 
very upper layers, eclipse measurements probe the atmosphere tangentially 
right at the cloud tops and thus enable observations with several times 
the integrated column density available in normal incidence observations. 
Thus the eclipse technique has the potential to reveal new details of 
the Jovian atmosphere. In order to obtain narrow band (40-360 K) eclipse 
measurements, a large aperture telescope and a fast response photo- 
electric spectrometer are necessary. This is true because the most 
informative portion of the eclipse light curve is that between 10th and 
15th magnitudes. In addition, this portion of the light curve is very 
short lived. 

The primary purpose of the Jovian Satellite Eclipse Study I was 
to obtain and document observations of five Galilean satellite eclipses 
in 1971 and to perform reduction of this data. The results of this 
study are contained in the present volume. 

Section II illustrates the effect of the various physical parameters 
on the eclipse light curve. In the next two sections the equipment and 
measurement procedures are detailed. The satellite absolute flux cali- 
brations follow. Details on the techniques employed in reducing the 
raw data to eclipse light curves, and the method for calculating theo- 
retical light curves are then presented. Finally the theoretical and 
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observational studies are combined and a number of preliminary con- 
clusions derived. A radius for Europa, and a Jovian atmosphere scale 
height are deduced. An aerosol extinction factor is derived from 
measurements at two wavelengths and the altitude of the cloud top is 
determined. The observed absorption features are tabulated and a 
number of additional data reduction tasks are outlined. 
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II. SATELLITE ECLIPSE DESCRIPTION 

In this section the significant physical parameters determining the 
shape of a Jovian satellite eclipse light curve are briefly introduced. 

For a detailed description, see Price (1970). From geometrical considera- 
tions, the best configuration for observing a Galilean satellite eclipse 
is quadrature. The phase angle (planetocentric angle between the Sun 
and Earth) is at its maximum, approximately 11° and the satellites enter 
and exit the shadow cone at the greatest angular distance from the limb 
of Jupiter. Ingress eclipse light curves are experimentally easier to 
observe in that the satellite may be acquired before the onset of the 
eclipse. Egress eclipse light curves are of interest because they enable 
study of the sunrise terminator but are more difficult to obtain because 
the satellite is not visible initially and off-set finding procedures 
must be used for acquisition of the satellite. Hence eclipses at western 
quadrature (before opposition) have been observed more frequently. Table 
II-l lists several physical data describing the satellites JI through JIV. 


TABLE II-l PHYSICAL DATA OF THE JOVIAN SATELLITES AND JUPITER* 


Name 

Magni tude 
At 

Opposi tion 
(%) 

Distance 
From 
Primary 
( Jovi an Radi i ) 

Si dereal 
Peri od 
(days) 

Radi us 
(km) 

lo (JI) 

+5.52 

5.92 

1 .769 

1670 

Europa (JII) 

+5.67 

9.40 

3.551 

1460 

Ganymede (Jill) 

+5.18 

15.00 

7.154 

2550 

Callisto (JIV) 

+6.20 

26.39 

16.69 

2360 

Jupi ter 
*Allen (1962) 

-2.50 

5.203 

(Au) 

11.862 

(years) 

71 ,350 
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During an eclipse a satellite first passes into the penumbra. The 
shape of the eclipse light curve during this portion of the eclipse is 
determined almost totally by the geometry of the problem i.e., by the 
sun and satellite diameters and distances from Jupiter. In the umbra 
the satellite would be totally dark if Jupiter had no atmosphere. (See 
cover picture.) However, Jupiter has an atmosphere which refracts the 
light so that the satellite remains visible after passing totallv into 
the umbra. The shape of this "refractive tail" on the eclipse light 
curve contains a great deal of information about the Jovian atmosphere. 
The remainder of this section is a qualitative description of the effect 
of each of the significant parameters on the light curves. Graphical 

data are employed. They were obtained in part from our own calculations 
performed in the manner described in Section VII and in part from Price 
(1970). 



Figure II-l Effect of Satellite 
Diameter on the Eclipse Light 
Curve (Price 1970) 


1) Satellite Radius. Figure II-l shows the effect of the satellite 

radius on the eclipse light curve. In the case of a point satellite, 
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the slope is steep but finite due to the finite diameter of the sun. 



TIME 


Figure II-2 Effect of Rayleigh 
Scattering for Various Wavelenaths 
(Price 1970) 


2) Rayleigh Scattering. One major extinction mechanism attenuating a 
ray as it passes tangentially through the Jovian atmosphere is 
Rayleigh scattering. Given the Hydrogen and Helium number densities 
versus altitude, the Rayleigh scattering cross sections, versus 
wavelength, are easily calculated. Figure 1 1-2 shows the effect of 
different amounts of Rayleigh scattering stemming from variation of 
the wavelength of the light in which the satellite is viewed. 

3) Scale Height. Figure II-3 shows the effect of the Jovian atmosphere 

scale height on the refractive tail. The angle of refraction through 
which a ray is bent depends inversely on the square root of the scale 
height, and on the closest approach elevation. Hence the shape of 
the refractive tail on the light curve is readily understood to be 
influenced by the scale height. 
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Figure II-4 Effect of Additional 
Aerosol Extinction on the Expected 
Refractive Tail 


Figure II-3 Effect of Scale Fleight 
on the Eclipse Light Curves for a 
Fixed Wavelength (Price 1970) 



4) Additional Extinction Due to Aerosols. It has been found necessary to 
invoke an additional extinction component to match our observed eclipse 
light curves with theoretical curves. If the ratio of aerosol and 
gaseous particles is constant throughout the atmosphere the effect of 
aerosol extinction has exactly the same form as that of Rayleigh 
scattering. See Figure II-4. The wavelength variation of this additional 
extinction component will yield the size parameter (ttx Particle 
diameter/X) variation of the aerosol particle extinction efficiency 
factors . 
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5) Additional Extinction Due To Molecular Absorption Bands. In particular 
bands additional extinction occurs due to absorption. The difference 
between the Rayleigh-plus-aerosol-continuum and the actual eclipse light 
curve in an absorption band enables determination of the absorbing 
specie number density, presuming the curve-of-growth is known. Figure 
II-5 illustrates this effect. 



Figure 1 1 -6 Effect of Cloud Top 
Cutoff on the Refractive Tail 


(6) Cloud Top Cutoff. A ray traversing tangentially through the atmosphere 
of Jupiter may get to a low enough altitude, at closest approach, to 
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be totally extinguished by the cloud tops. This terminates the refractive 
tail and is illustrated in Figure II-6. 

In principal, the ray might also be terminated by critical refraction-- 
the situation in which the ray is refracted with a smaller radius of 
curvature than the planet radius. The observational effect would be 
exactly the same as with cloud top cutoff. 



TIME 

Figure II-7 Effect of the Figure 1 1 -8 The Shadow Chord 

Albedo Change Producing Length Determined from Ingress/ 

Anomolous Post Reappearance Egress Observations 

B righting 

7) Satellite Albedo and Anomalous Reappearance Brightening. If the satellite 
reflectivity should increase during immersion in the shadow the satellite 
might reappear with a greater brightness than prior to disappearance. 

In principal, detection of such brightening is possible from reappearance 
data. See Figure II-7. 

8) Jovian Figure Determination. Figure II-8 illustrates how knowledge of 
the precise duration of a satellite eclipse may lead to a determination 
of the chord length of the shadow cone, and hence of the Jovian disk. 
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III. EQUIPMENT 

The 200-inch Hale telescope and the 33-channel photoelectric spec- 
trometer were used in this study to measure the eclipse light curves. 

The spectrometer optical system and associated electronic equipment 
are described by Oke (1969). The spectrometer is designed for operation 
at the Cassegrain focus. The data were recorded on punched cards and 
printed paper tape. Selected wavelength counters were also recorded on 
video tape. The punched card mode was used during the non time-critical 
portions of the measurements. The printed paper tape mode was employed 
at all times, and the video tape mode was employed only during the 
actual eclipses. An audio tape machine recorded voice remarks from the 
data room and from the observer in the Cassegrain cage on one channel 
and WWV on the other channel. Audio recordings were made at all times, 
including set-up and calibration periods. 

A. Multichannel Photoelectric Spectrometer 

Only a brief description of the spectrometer will be given here. 

o 

The spectrometer operates in the wavelength region 3100 to 11,000 A. 
Although there are 33 photomultiplier channels available, only 20 were 
used because only 20 sets of dual pulse counters were available. The 
wavelength interval is divided into two regions, 3100 to 5600 K (ref- 
fered to as blue)and 5600 to 11,000 K (referred to as red). In the 

o 

blue region, the 17 photomultipliers have acceptance bands 160 A wide; 
while in the red region, the 16 photomultipliers have acceptance bands 
360 K wide. In the focal plane of the spectrometer camera, a mask or 
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decker can be inserted so that any restricted wavelength interval 
within each of the bands in the red or blue can be measured. Various 
standard "Oke-decker" patterns are available in the red and blue with 
20, 40, 80, 160 and 360 A slots. The grating can also be rotated so 
that a spectral scan is obtained. For the eclipse measurements, 
special red and blue "Boeing-deckers" were constructed. On February 8, 
the specially designed Boeing red decker was used but the blue decker 

o 

could not be inserted into the spectrometer and the 160 A band, blue 
Oke decker #4 was used during that eclipse. For the March and April 
eclipses the Boeing deckers were used in both the red and blue. 

There are dual entrance apertures to the spectrometer separted by 
40 arc seconds. A choice of entrance aperture diameters is available 
as shown in Table III-l. A chopper behind these entrance apertures 
alternately admits light from first one and then the other aperture 
at 30 Hz (see Figure VI-2), and the counts from each aperture are 
electronically separated. Various integration times may be selected. 

In addition to the counts, which are printed out for each aperture and 
the 20 selected channels, other information such as integration time, 
the grating position, entrance aperture number, decker number, right 
ascension, declination, date , sidereal time, civil time, and certain 
alphameric identifications can be selected for printing and/or punching. 

Controls for the multichannel spectrometer and the data systems 
are located both in the data room and at the Cassegrain focus. Nixie 
tube displays are located in the data room for the 20 selected channels 
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TABLE III-l SPECTROMETER INSTRUMENT ENTRANCE APERTURES 


NO. 

DIAMETER 
(Sec arc) 

DIAMETER 

(mm) 

0 

1 

2.03 

0.79 

2 

3.58 

1.39 

3 

5.04 

1 .96 

4 

6.94 

2.70 

5 

9.92 

3.86 

6 

14.05 

5.46 

7 

20.00 

7.77 


Separation of holes = 40.00 sec arc = 15.57 mm 

Except for No.l, apertures increase in area by 
very nearly factors of 2.0 
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along with time, telescope position, spectrometer and data system status. 
An abbreviated display is located at the Cassegrain focus. 

B. Additional -Recording Systems 

A two channel audio tape recorder was employed to keep a log of all 
checkout and calibration procedures which began about midnight for each 
eclipse. One channel recorded the output of the intercommunication 
system between the Cassegrain cage and the data room, plus all the ac- 
tivities in the data room. A second channel was used to record WWV 

received at 5 Me. This system was found particularly useful in that 

a few deviations were made from the planned observing sequence of 
events during the experiment which were not written down at the time. 

In addition to audio recording, a TV camera and TV recording system 
was used. The most important channels in the red and blue, and civil 
and sideral time were displayed on Nixie tubes within the field of 
view of the TV system. During the actual eclipses, only the printed 
paper tape data recording mode was used because the card punch was 
slower, requiring ~15 seconds to punch a series of observations 
after each 10 second integration period. The printed paper tape 
required only 8 seconds to print after each 10 second integration 

period. In the event of a printer system failure the video tape backup 

data system would have enabled data retrieval in several important 
channels. In addition, a standby printer was available. During the 3 
nights of observation when the tape printer was in continual operation 
there was no failure of this or any other electronic system. 
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IV. MEASUREMENTS 


Introduction 

Four satellite ingress and one satellite egress eclipses were 
observed at 20 different wavelength bands making a total of 100 eclipse 
light curves. The 20 wavelength bands were slightly different for the 
three observation periods. The choice of continuum bands and exoected 
features were, in general, determined from the literature (Kiess et al . , 
I960; Spinrad, 1963; Ov/en , 1967; Spinrad and Trafton, 1963; and Taylor, 
1965). 

After the discriminator curves of the photomultiplier tubes were 
checked for each channel, an alternate choice of channels was sometimes 
necessary because of a bad tube. Several spare wavelength band passes 
were built into the deckers in view of this possibility. Table IV-1 
lists the satellite eclipses that were observed in 1971. Each observa- 
tion included observations of the star HD 140283. This star was employed 
as a standard star (Oke, 1971). Typically a measurement sequence con- 
sisted of the following events: reduce the telescope aperture to nomi- 

nally 100-inches, next observe the calibration star, then observe the 
eclipse until the light level reaches a predetermined value, then open 
the telescope to full aperture (200-inches) and observe the remainder of 
the eclipse. Finally the standard star was again observed at the 200 
inch aperture. 

The dust cover mechanism on the 200-inch telescope can be used to 
reduce the collecting area of the primary mirror. It was necessarv to 
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TABLE IV-1 SATELLITE ECLIPSE OBSERVATIONS OF JI, JII, & Jill 


SPECTROMETER 
. ENTRANCE 
AIR"^ APERTURE 

DATE (1971) SATELLITE ECLIPSE TIME (PST) Rx R * MASS (ARC SEC) 


February 8 

JII (Europa) 

Ingress 

5:39 

am 

-2.6 

-0.5 

1.8 

6.94 

March 10 

Jill (Ganymede) 

Ingress 

2:56 

am 

-3.2 

-0.8 

2.0 

14.05 



Egress 

5:10 

am 

-2.0 

-0.8 

1.7 

14.05 

April 6 

JII (Europa) 

Ingress 

2:22 

am 

-2.1 

-0.5 

1.8 

5.04 


JI (lo) 

Ingress 

4:02 

am 

-1.8 

-0.3 

1.8 

3.58 


* 

Rx = Distance in Jovian radii from disk center. Eastward (parallel to equator of 
Jupiter) to eclipse location 

Ry = Distance in Jovian radii from disk center. Northward (parallel to axis of 
Jupiter) to eclipse location 

^ Air Mass = Air mass of Jupiter from Palomar at eclipse time 
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reduce the telescope aperture to a value between 90 and 100 inches 
during the early stages of the eclipse (m^ ~ 5) to prevent satura- 
tion of the photomultiplier tubes in the multichannel spectrometer. 

When the brightness of the eclipsing satellite was sufficiently 
reduced, the dust cover was opened completely permitting measurements 
of the satellite to continue to much lower light levels (m^ ~ 15). 

The brightness of HD 140283 is m^ = 7.3 so that observation at full 
aperture presented no saturation problem. 

Calibrations using the standard star were performed each night 
using the Oke red #2 and blue #1 deckers and the Boeing red and blue 
deckers. Spectral scans were made once each night in 10 Angstrom 
steps in the blue and 20 in the red. Subsequent calibrations during 
the night were then made with the grating fixed in the nominal posi- 
tion (500). 

For each series of observations a judgment of the seeing conditions 
was made and the smallest usable spectrometer instrument aperture selected 
(see Table III-l). The instrument aperture setting was reduced when 
seeing conditions permitted. 

In order to evaluate the sky background around the satellite, 
periodically during the eclipse, the object aperture was: 1) moved 

off the satellite toward Jupiter, 2) moved back on the satellite. 
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3) moved off on the side away from Jupiter, 4) then moved back on. 
Measurements were made at all four positions. Extrapolation between 
the two readings from the object-aperture adjacent to the satellite 
could then give the background sky at the satellite position. This 
procedure is referred to as a "spatial scan" in subsequent descrip- 
tions . 

A. February 8, 1971 Observation of Europa 

This was the first eclipse observation so that particular care 
was taken in calibration procedures and experimental techniques. Table 
IV-2 lists the sequence of events. Here a few explanations of the 
abbreviations are given: the oke decker Red #2 has 40 Angstrom slits 

for each channel. The oke decker Blue #4 has 160 Angstrom slits (wide 
open) for each channel. The term R/s aperture #s refers to the orien- 
tation of the two spectrometer instrument apertures. The number refers 
to one of the apertures listed in Table III-l. The object-aperture is 
the one in which the object is being observed. For example, in Table 
IV-2, at time 00:57, Observe Uranus^ grating (400,20,600) S Aperture 
#5, means: with the telescope closed down to 100 inches and Oke decker 

Red #2 and Blue #1 in place, measurements at 10 second integration 
intervals were taken with Uranus in the south aperture. The instrument 
aperture diameter was 9.92 arc seconds. The grating was set at reading 
400 (100 Angstroms below the nominal position) and a data count taken. 
The grating was then rotated to reading 420 and a data count taken. 
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TABLE IV-2 SEQUENCE OF EVENTS FOR FEBRUARY 8, 1971 OBSERVATIONS 

OF JII (EUROPA) 


TIME 

00:00 

00:57 

01 :25 

01 :36 
02:04 

02:24 

02:54 

03:02 


LI EVENT 

Check out equipment, final selection of channels 
• Close telescope to 100-inch (1st setting) 

Set on Uranus, estimate seeing, determine spectrometer 
entrance aperture 
Obtain 16 mm movie sequence 
Oke decker Red #2, Blue #1, integration time 10 sec 

Observe Uranus , grating (400,20,600)* S aperture #5 
Observe Uranus , grating (600,20,400), N aperture #5 
Boeing decker Red #0, Oke decker Blue #4 

Observe Uranus , grating (480,10,520) N aperture #5 

Observe Uranus , grating (520,10,480) S aperture #5 

Rehearse satellite eclipse sequence with printer 
only, video tape recorder, etc. 

Oke decker Red #2, Blue #1 

Observe Uranus , grating (400,20,600) S aperture #5 

Observe Uranus , grating (600,20,400) N aperture #5 

Boeing decker Red #0, Oke decker Blue #4 
Practice satellite eclipse sequence 
Observe HD 140283 , grating (480,10,520) N aperture #5 
Observe HD 140283 , grating (520,10,480) S aperture #5 
Oke decker Red #2, Blue #1 
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TIME (PST) 


EVENTS 


Observe HD 140283 , grating (400, 10,600) N aperture #5 

Observe HD 140283 , grating (600, 10,400) S aperture #5 

03:26 • Open telescope to 200-inch 

Boeing decker Red #0, Oke decker Blue #4 

Observe HD 140283 , grating (500) N/S aperture #5 
04:16 •Close telescope to 100-inch (2nd setting) 

Observe HD 140283 , grating (500) N aperture #5 

Observe HD 140283 , grating (500) S aperture #5 

04:19 Set on JII (Europe) ^ estimate cage rotation, move to 

zenith rotate cage, 16 mm sequence, Boeing decker 
still in position, grating set at (500) 

Observe JII, practice eclipse routine, on JII; off 
JII, toward Jupiter; on JII; off JII away from Jupiter. 
04:30 Observe HD 140283 , grating (500) S aperture #4 

Observe HD 140283 grating (500) N aperture #4 
04:40 Set on JII, practice eclipse sequence 

05:01 Observer HD 140283 , grating (500) S/N aperture #4 

05:10 Set on JII, Boeing decker, grating (500), 100-inch 

(2nd setting) 

JII Eclipse Sequence 10 sec integration time, S 
aperture #4 

Repeat spatial scan-observe sequence: 

(a) 0*^ 0^ start a single counting sequence 
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TIME (PST) 


EVENT 


m Q 

(b) 0 18 move toward Jupiter one aperture 

diameter (14 arc second) 

(c) O"’ 28^ count 

(d) 46^ move back on JII 

(e) 0*^ 56^ count 

(f) 14^ move off JII away from Jupiter one aperture 
diameter (14" arc second) 

(g) 1*^ 24^ count 

(h) 1*^ 42^ move back on JII 

(i) l"’ 52^ count 

(o) 2*^ 10^ repeat count approximately every 18 seconds 

m s 

(k) 6 56 repeat sequence (a) through (j) 

05:30 Repeat above sequence until counts on channel 17 

drop below 250,000 and then, 

05:40 Open telescope to 200-inch aperture 

Continue to observe JII 
05:53 End eclipse observation 

05:54 Observe HD 140283 grating (500) S/N aperture #4 

05:58 Conclude measurements 

if 

(grating 400,20,600) = Move grating in steps of 20 from 400 to 600. 

Some notation used for steos of 10 between other limits. 

ifi( 

For telescope apertures other than full (200-inch) it is not 
possible to set with precision. Calibrations must be performed 
for each aperture less than 200. The notation 1st, 2nd . . . 
setting is employed in all Tables. 
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This procedure was continued up to grating setting 600. For example, 
on Channel #5 at grating reading 460 the wavelength interval measured 
was XX 3660 to 3680 Angstroms. 

Spectral scans were made of Uranus using the two sets of deckers. 
The purpose was to observe the H2 dipole line and check the wavelength 
position of the Boeing deckers. 

The wavelength bands observed February 8 are listed in Table IV-3. 
Several column headings are explained: 

CHAN. NO. is the photomultiplier channel sequence number from 
shortest wavelength (blue) to longest wavelength (red). 

SEQUE. NO. is the sequence number assigned to the channel. 

DECKER USED FOR OBSERV. is the decker type used for the obser- 
vation. The number is the physical order in which the Nixie 
tube display was arranged in the data room and the order in 
which that data were printed out on the printed paper tape. 

B. March 10, 1971 Observation of Ganymede 

This was the second eclipse observation so that experience from 
the earlier observation resulted in some changes in experimental pro- 
cedures. Table IV-4 lists the sequence of events. In addition to a 
spectral scan of Uranus, one of Neptune was also accomplished to further 
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TABLE IV-4 SEQUENCE OF EVENTS FOR MARCH 10, 1971 GALILLEAN 
SATELLITE ECLIPSE OBSERVATIONS Jill (GANYMEDE) 

TIME (PST) EVENT 

00:00 Check equipment, ice spectrometer, rotate spectrom- 

eter, dark current readings 
00:30 • Close telescope to 100-inch (1st setting) 

Oke decker Red #2, Blue #1, integration time 10 sec 
00:43 Observe Uranus , grating (400,20,600) N/S aperture #6 

00:55 Observe Sl^ 10 arc second N of Uranus, grating (500) 

01:00 Boeing decker Red #0, Blue #0 

Observe Uranus , grating (460,20,500) N/S aperture #6 
Observe Sky , 10 arc second N of Uranus, grating (500) 
01:28 Observe HD 140283 , grating (500) N/S aperture #6 

Oke decker Red #2, Blue #1 

Observe HD 140283 , grating (400,20,600) N/S aperture #6 

• Open telescope to 200-inch 

01:51 Observe HD 140283 , grating (500), N/S aperture #6 

Boeing decker Red #0, Blue #0 

01:55 Observe HD 140283 , grating (500), N/S aperture #6 

• Close telescope 100-inch (2nd setting) 

02:00 Observe HD 140283 , grating (500), N/S aperture #6 

Check rotation of spectrometer on Jill 
02:11 Observe Jill grating (500), N/S aperture #6 
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TIME 

02:13 

02:17 

02:24 

02:28 


03:02 


03:13 


03:35 


03:38 

03:41 


LI EVENT 

•Close telescope to 90-inch (1st setting) 

Observe Jill grating (500), N/S aperture #6 
Observe HD 140283 , grating (500), N/S aperture #6 
Observe Jill , grating (500) 

Jill Eclipse Sequence (background determination included) 

•Every 5 min move off Jill 14 arc second toward Jupiter, 
then count, move back to Jill, 

then count, move off Jill 14 arc second away from Jupiter, 
then count, move back to Jill, 
then count, continue to count 
Note when count down by a factor of 10 in channels 22 & 11 
•Open telescope to 200-inch 

Continue counting well into umbra 
Repeat background determination 

Oke decker Red #2, Blue #1 

Observe Jill (not now visible) grating (400,20,600) 

S aperture #6 

Boeing decker Red #0, Blue #0 

Observe scattered light , grating (500) aperture #6 
Move 20 arc second from Jupiter's limb, count 
Move in 20 arc second steps away from Jupiter, count 
repeat 4 times 

Observe HD 140283 , grating (500), N/S aperture #6 
Close telescope to 100-inch (3rd setting) 
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TIME 

03:50 

03:58 

04:27 

04:33 

04:37 

04:39 

04:49 

05:12 

05:15 

05:18 

05:21 

05:23 

05:24 

05:26 

05:54 


n EVENT 

Observe HD 140283 , grating (500), N/S aperture #6 
Observe Neptune , grating (460,20,540), N/S aperture #6 
Oke decker Red #2, Blue #1 

Observe Neptune , grating (460,20,540), N/S aperture #6 
Observe HD 140283 , grating (500), N/S aperture #6 
Boeing decker Red #0, Blue #0 

Observe HD 140283 , grating (500), N/S aperture #6 
• Open telescope to 200-Inch 

Observe HD 140283 , grating (500), N/S aperture #6 
•Close telescope to 100-Inch (4th setting) 

Observe HD 140283 , grating (500), N/S aperture #6 
Observe Jill , grating (500), S aperture #6 

Jill Eclipse Egress Sequence 

Every 5 m1n (same as 02:30 sequence) 

Observe Jill , grating (500), S aperture #5, count once 
Observe HD 140283 , grating (500), S/N aperture #6 
•Close telescope to 90-inch (2nd setting) 

Observe HD 140283 , grating (500), S aperture #5, count 
twice 

Observe Jill , grating (500), S aperture #5 
Conclude Jill Egress sequence 

Observe Scattered light radially away from Jupiter 
Take dark count, rerun discriminator curves 
Conclude measurements 
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verify the location of the H 2 dipole band. Since Ganymede (Jill) was 
brighter than JII the telescope aperture had to be reduced to 90 inches 
for the initial portion of the eclipse. During the totality phase of 
the satellite eclipse a search was made for any extraneous emission 
from Jupiter's nighttime side by scanning Jill in wavelength. Several 
measurements were conducted to evaluate the effects of scattered light 
from Jupiter. 

The telescope was set at the 100-inch aperture setting to observe 
the egress of Ganymede instead of starting at full aperture (200-inch). 
This was done to avoid loss of time in closing the telescope down to a 
smaller aperture when the satellite brightness reached a level that 
would saturate the photomultipliers. Also, because of twilight, a 
calibration of the standard star might have been difficult. An analysis 
of the data has shown observation well into twilight is possible so that 
subsequent observations can be made for longer periods. 

Table IV-5 lists the wavelength band used for Ganymede. The 
Boeing decker Blue #0 was altered to fit the spectrometer. For this 
observation, channel number 16 was 60 Angstroms wide to observe the 
NH 3 band and channel number 17 was 120 Angstroms wide to better fit 
the continuum than the channel used in February. Channel number 18, 
which had not been functioning previously was used so that the con- 
tinuum and CH 4 could be observed together in a 360 Angstrom band. 

This was the only night that channel number 18 was used. The order 
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of the physical display of channels was changed so that channels 
number 3, 8, 11, 18, 20, 22, 31 would be in the field of view of the 
TV recorder system. This was done so that in the event of a printer 
failure the most significant data could still be obtained. 

C. April 6, 1971 Observation of lo and Europa 

April 6 provided the unique opportunity to observe two eclipses, 
an ingress of JI and an ingress of JII. Table IV-6 lists the sequence 
of events for the two eclipses. Additional measurements were made to 
evaluate the scattered light effect from Jupiter and the effect of 
different spectrometer instrument apertures. 

At the conclusion of the second eclipse observation, a series of 
twilight measurements were made at -23° declination. The purpose of 
these measurements was to determine the sky brightness under the 
same conditions that will prevail during an egress of JIV in 1972. 

The wavelength bands used for the April 6 observation are listed 
in Table IV-7. Several changes were made from the February and March 
wavelength bands. Channel number 13, a continuum plus CH^ region^ was 
used and the bluest band, 3160 to 3320 Angstroms, was eliminated. 
Channel number 18, 5640 to 6000 Angstroms, was also eliminated. 

In summary, all the planned observations were made. Seeing con- 
ditions were not always ideal and are discussed in Section VI, Data 
Reduction. There was no malfunction of any electronic equipment. 
During the twilight measurements the wind screen was found to be 
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TABLE IV-6 SEQUENCE OF EVENTS FOR APRIL 6, 1971 
GALILLEAN SATELLITE ECLIPSE OBSERVATIONS JI (10) AND JII (EUROPA) 

TIME (PST) EVENT 


00:00 

00:12 


01 :10 
01 :32 


01 :39 


01:48 
01:53 
01 :59 


Check equipment, ice spectrometer, rotate spectrom- 
eter, dark current readings 

• Open telescope to 200-inch 

Oke decker Red #2, Blue #1 

Observe HD 140283 , grating (370,20,600) S/N aperture #3 
Observe HD 156074, grating (370,390,420,500,550,570,580, 
600) S/N aperture #3 
Boeing decker Red #0, Blue #0 

Observe Sky near JII, between JI and JIV, aperture #3 
Check rotation 

• Close telescope to 90-inches (1st setting) 

Observe ^I_, grating (460, 20, 540) S aperture #2, 

#3, #1 

Observe JII , grating (500) S aperture #3 
Observe HD 140283 , grating (500) S/N aperture #3 
JII Eclipse Sequence , grating (500) S aperture #3 

Every 5 min move off JII toward Jupiter, count 
move back on JII count 
move off JII away from Jupiter, count 
move back on JII, count 
continue count 
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TIME (PST) 


EVENT 


Note when count down by a factor of 8 in channels 
22 and 11 

02:22 * 0pen telescope to 200-inch 

Continue counting into umbra 
Repeat background determination 
02:43 Complete JII eclipse measurement 

02:46 Observe HD 140283 , grating (500) N/S aperture #3 

03:02 Observe scattered light from Jupiter, grating (500) 

Set aperture #1 move to center Jupiter 

Move in RA one Jupiter radius off W limb, count 

Set aperture #2, count 

Move in RA two Jupiter radii off W limb, count 
Move in RA three Jupiter radii off W limb, count 
Set aperture #3, count 

Move in RA four Jupiter radii off W limb, count 
Set aperture #4, count 

Move in RA five Jupiter radii off W limb, count 
Move in RA six Jupiter radii off W limb, count 
03:20 Move in RA three Jupiter radii off W limb 

aperture #1 , count 
aperture #2, count 
aperture #3, count 
aperture #4, count 
aperture #5, count 
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TIME (PST) 


EVENT 


Set aperture #6, count 

Set aperture #7, count 

03:24 Dark count, 50 second integration time 

03:27 Sky evalution 

Aperture #2, count 
Aperture #3, count 
Aperture #4, count 
Aperture #6, count 
Aperture #7, count 

Long Sky count near HD 140283, grating (500) 
aperture #3, 100 sec integration 
03:36 • Close telescope to 90-inch (2nd setting) 

Observe HD 140283 , grating (500) N/S aperture #3 

03:45 J1 Eclipse Sequence , grating (500) S aperture #2 

Every 5 min repeat off and back sequence (see 01:59 
position for details) 

Note when count down by factor of 4 in channels 22 
and 11, then 

04:01 • Open telescope to 200-inch 

Continue counting into umbra 
04:16 Complete JI eclipse measurements 

Evaluate Sky near HD 140283 

Sky , grating (500) aperture #7, 100 sec integration 
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TIME (PST) 


EVENT 


04:25 

Observe HD 140283, grating (500) N/S aperture #2 

04:31 

Evaluate twilight background at -23° declination, 
grating (500) N/S aperture #7, continue to count 

04:55 

Continue count 

04:58 

Aperture #6, continue count 

05:01 

Aperture #5, continue count 

05:04 

Aperture #4, continue count 

05:11 

Aperture #2, continue count 

05:12 

Aperture #4, continue count 

05:17 , 

Aperture #1, continue count 

05:25 

Completed measurement, wind screen found to be in 
field of view at this time 

06:36 

Concluded measurements 
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partially blocking the field of view and the observation was concluded. 
However, the eclipse data sought were apparently obtained satisfac- 
torily. So far, we have not been able to ascertain the precise time 
of wind screen malfunction from the audio tapes. It is felt that it 
failed after viewing the standard star at the conclusion of the lo 
observations . 
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V. CALIBRATIONS 


The initial point on each of the light curves obtained during 
three eclipses has been reduced to an absolute monochromatic flux 
above the earth's atmosphere. The three eclipses so analyzed were 
those of lo, April 6, 1971; Europe, April 6, 1971; and the Ganymede 
ingress, March 10,1971. The reductions were performed with a code 
kindly made available to the authors by Dr. J. B. Oke. The results 
are presented in Table V-1 , 2, and 3. ' It is estimated that the 
standard deviation of the data in the Magnitude Column does not exceed 
+0.05 magnitude, except for the data marked with the colon (:), 
indicating more uncertainty. 
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TABLE V-1 


ABSOLUTE FLUX ABOVE THE EARTH'S ATMOSPHERE 
lo, April 6, 1971 


BANDPASS CENTER MONOCHROMATIC FLUX, LOG MONOCHROMATIC FLUX 

^ MAGNITUDES . . , 

(A) (-2.5 log fy -48.60) (^og-iQ erg cm"^sec" Hz" ) 


3400 

10.60 : 

-23.68 

3720 

10.03 : 

-23.45 

4200 

7.40 

-22.40 

4680 

6.29 

-21.96 

5000 

5.88 ; 

-21.79 

5320 

5.70 

-21.72 

5530 

5.75 

-21.74 

5660 

5.54 

-21.65 

6325 

5.37 

-21.59 

6420 

5.30 

-21.56 

6850 

4.71 

-21.32 

7265 

5.02 

-21.45 

7495 

4.95 : 

-21.42 

7978 

5.23 : 

-21.53 

8260 

5.11 

-21 .48 

9115 

5.12 

-21.49 

9260 

5.32 

-21.57 

10140 

5.05 

-21.46 

10500 

5.06 

-21.46 

10860 

5.11 

-21.48 


36 



TABLE V-2 

ABSOLUTE FLUX ABOVE THE EARTH'S ATMOSPHERE 
Europa, April 6, 1971 


BANDPASS CENTER MONOCHROMATIC FLUX LOG MONOCHROMATIC FLUX 

MAGNITUDES 

(A) (-2.5 log -48.60) (log^o erg cm"^sec"^Hz“^ ) 


3400 

8.65 

3720 

7.77 

4200 

6.53 

4680 

6.05 

5000 

5.83 

5320 

5.71 

5530 

5.76 

5660 

5.55 

6325 

5.49 

6420 

5.43 

6850 

4.96 

7265 

5.21 

7495 

5.43 

7978 

5.15 

8260 

5.29 

9115 

5.31 

9260 

5.51 

10140 

5.24 

10500 

5.30 

10860 

5.33 


-22.90 : 
-22.55 : 
-22.05 
- 21.86 
-21.77 : 
-21.72 
-21.74 
- 21.66 
-21.64 
-21.61 
-21.42 
-21.52 
-21.61 : 
-21.50 : 
-21.55 
-21.56 
-21.64 
-21.53 
-21.56 
-21.57 
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TABLE V-3 

ABSOLUTE FLUX ABOVE THE EARTH'S ATMOSPHERE 
Ganymede, March 10, 1971 

BANDPASS CENTER MONOCHROMATIC FLUX LOG MONOCHROMATIC FLUX 

MAGNITUDES 


(A) 

(-2.5 log -48.60) 

{log-jQ f^, erg cm"^sec“^Hz 

3240 

7.84 : 

-22.58 r 

3400 

7.52 : 

-22.45 : 

3720 

6.98 ; 

-22.23 : 

4200 

5.94 

-21.82 

4680 

5.51 

-21.64 

5320 

5.16 

-21.50 

5530 

5.22 

-21.53 

5660 

5.02 

-21 .45 

5820 

5.00 

-21.44 

6420 

4.91 

-21.40 

6850 

4.42 

-21.21 

7265 

4.71 

-21.32 

7495 

4.93 : 

-21.41 : 

7978 

4.65 : 

-21.30 : 

8260 

4.79 

-21.35 

9115 

4.85 

-21.38 

9260 

5.11 

-21.48 

10140 

4.76 

-21.34 

10500 

4.81 

-21.36 

10860 

4.88 

-21.39 
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VI. DATA REDUCTION 


Introducti on 

In this section the technique used to obtain the satellite eclipse 
light curves from the raw data is described. A detailed description of 
the measurement sequence was presented in Section IV. Briefly, the 
telescope was stopped down to nominally 100 inches and after observing 
the standard star, observation of the satellite was begun. Prior to 
initiation of the eclipse ingress, occasional sky data were obtained 
just east and just west of the satellite position. After eclipse 
ingress had proceeded to a predetermined count level the telescope was 
opened to the full 200-inch aperture. Observations were continued well 
after the satellite was indistinguishable against the background sky. 

The standard star was again observed. A series of measurements was 
obtained approximately every 18 seconds consisting of sky background 
counts and satellite plus background counts in 20 wavelength bands. 

These data were then reduced to obtain the eclipse light curves. 

An eclipse light curve is a plot of the darkening, in magnitudes, relative 
to the pre-eclipse brightness, versus universal time (see Figure Vl-la). 
The first step in the reduction procedure was removal of the background 
sky brightness from the object (sky background plus satellite) data. 

Two methods have been employed. One is to determine a "balance factor", 
the relative sensitivity of the two spectrometer apertures. See Figure 
VI-2. Using this balance factor, the count obtained from the object 
aperture was adjusted by subtracting a corrected sky brightness. The 
residual is the corrected object count, i.e., satellite brightness. 
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TIME ORIGIN, MARCH 10. 1971 


The drawback to this dual aperture sky removal technique is that the 
residual true object count obtained suffers from the fluctuations in 
both the object and sky counts relative to their mean, and no advantage 
is taken of the possibility to average a number of sky counts and hence 
to obtain an improved corrected object brightness. The lack of these 
refinements is of no consequence during the brighter portions of the 
eclipse, but during the refractive tail portion it is desirable to utilize 
any refinement possible. Hence the other method employed, the curve fit 
sky removal technique , utilizes a curve, least squares fitted through 
the object counts observed after the satellite had become invisible to 
the spectrometer. This curve was then extrapolated to earlier times 
prior to satellite disappearance. By subtracting this extrapolated sky 
count from the object count the corrected object count was obtained as 
illustrated in Figure Vl-lb. This method was employed to correct the 
important refractive tail portion of the light curves obtained during 
times of good sky stability. 

The general subject of sky removal from the data has introduced 
the necessity of determining balance factors. Hence this subject is 
discussed below in more detail. There is also the necessity of determining 
the change in effective sensitivity of the system introduced by using 
different telescope collecting areas. This subject is discussed below. 

Next there is a summary of the time scale correction procedure. A 
discussion of errors follows. Finally, all the elements of the data 
reduction procedure are summarized. 
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1. Balance Factors 


The primary advantage of a dual aperture spectrometer is that simultaneous 
observations may be made both of an object and of the adjacent sky background. 
The counts received in the object aperture also include counts from the 
sky background. In principle, the sky background may be removed by 
subtraction. However, there is always the possibility of an imbalance 
or difference in response when a signal of the same strength is observed 
first in one aperture and then in the other because of sensitivity 
variations over the photo tubes, geometrical effects, etc. In the case 
of the present experiment, there is an additional problem. There is a 
difference of sky background in the two apertures occasioned by the slight 
difference in scattered light which arises from the different distances, 
r.| and or the two apertures from the nearby, very bright limb of 
Jupiter (see Figure VI-2). Initially, r.j and were approximately the 
same, but they changed and at a different rate during the typically 
40-minute observation sequence. Therefore, a balance factor is required 
that is unique not only to the spectrometer wavelength channel but also 
to the particular observation. 

The procedures employed to obtain these balance factors and the 
corrected object counts were these. First an average was obtained for 
the object to sky ratio over 10 or more counts a little after the 
satellite was judged to be indistinguishable against the sky background -- 
"satellite-disappearance". Then the corrected sky count was subtracted 
from the object count to yield the true object count. The notation 
employed in these calculations is listed in Table VI-1. 
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Table IV-1 Notation Used for Balance Factor Calculations 


Subscript i - denotes the count, i=l,.,.,100+ depending on the number of 
integrations performed. 

Subscript j - denotes the wavelength band, j=l,...,20. 

Superscript k - denotes an iteration in calculating the corrected object 
count or balance factor. 

0— - denotes the number of counts observed in the object aperture 
on the ith observation in the jth wavelength channel 
S. . - simarly denotes the number of counts in the sky aperture 

* J 
1^ 

C. . - denotes the corrected object counts on the kth iteration 

* J 

BF. - denotes the balance factor for the jth wavelength channel 

J 

obtained on the kth iteration. 

CF^. - a correction factor, discussed below, for all channels on 
the ith count and obtained on the kth iteration. 

ML,NH - first and last values of i in the interval in which the 
average balance factor is obtained. 


Hence 



S. . X BF*^ 
ij J 


(VI-1) 


where 


BF*^ 

J 


NH 

E 


0. 


NL 


1J 


pk-1 

U ^ij 
(NH-NL+1) 


(VI-2) 
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Note that, in principle, the guess as to "satellite-disappearance" 

may have been early. Effort was made to make the guess as soon as possible 

after disappearance so that the different rates of change of r^ and 

would be insignificant. Hence, in determining the balance factor, two 

k *• 1 0 

iterations were used. On the first, C . . = C . . was set equal to 0. 

* J I J 

The first and second iteration balance factors differed by only a few 

5 

parts in 10 , typically. 

An additional phenomenon was observed in reducing the data. The 
object/sky ratio in each channel varied systematically relative to the 
channel average, from one observation to the next after the satellite 
had disappeared. This suggested that a further improvement in the 
correction procedure might be employed in which the sky reading, on any 
given observation, would be corrected by an additional factor determined 
from al 1 the wavelength channels. This correction factor, CF, essentially 
notes how the particular object/sky ratio on the ith observation and 
averaged over all 20 channels differed from the mean object/sky ratio 
obtained from averaging over all 20 channels and over the interval 
i=NL to NH. Quantitatively: 


20 

S ■ 

j = l 

(«ij - 

s.. 

jj V ij X ^ 

(VI-3) 

NH 

E 

i=NL 

^ 20 

0 . . - o':. ' V s. . 

1J 1J TJ 


S. (NH-NL+l) j=l 

' %J 



This last factor in VI-3 weights the individual channel 0. ./S. . 

1 J * J 

ratios according to the number of counts. The denominator in the first 

1/ 

factor on the right side of equation VI-3 normalizes CF^ to unity. 
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on the average. Equation VI-1, which gives the corrected object counts, 
then becomes 

. = 0. . - S. . X BF*^ X (VI-4) 

ij ij TJ J 1 

To summarize, a two-step iteration procedure was employed. The 
C.. 's were determined from equation VI-4 using balance factors obtained 
with equation VI-2 and correction factors obtained with equation VI-3. 

2. Telescope Aperture Area Ratios 

In order to relate the data obtained prior to opening the telescope to 
200-inch to the data obtained after, it was necessary to determine the 
effective system sensitivity change. This can be done by comparing the 
corrected count rate obtained when observing a standard star before and 
after fully opening the dust cover. This procedure was used in reducing 
the February ingress data. However, an alternate procedure was employed 
with the March and April data. The ratio was adjusted until the eclipse 
light curves manifest a smooth continuous fall off. This alternate 
technique was employed because the first technique did not give rise to 
smooth light curves in all cases. It is suspected that the spectrometer 
responds differently to an extended source, the satellite image, than 
to a point source, a standard star. It was also noted that the effective 
collecting area ratio is slightly channel (wave length) dependent. 

Careful inspection of the curves in Part Two for the April ingress of 
Europa evidences this effect. However, because the magnitude of the 
effect is significantly less than the uncertainty in the data points 
in the refractive tail, further correction of the area ratio factor for 
each channel was not done for this report. 
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3. Time Scale Conversion 

The time scale for the eclipse light curves was obtained in the following 
manner. Prior to beginning the data gathering sequence, a block of 
spectrometer and system configuration information was printed out (ID card). 
This data included both Pacific Standard Time and Sidereal Time. Each 
night, PST was set by WWV to 1/10 of a second. The sidereal time was 
also printed out with each set of count data. The time printed was the 
sidereal time at the conclusion of the 10.007 second integration period. 
These sidereal times were converted to UT using the UT equivalent of the 
initial PST and the UT interval equivalent of the sidereal time interval 
since the initial instrument configuration sidereal time. In this manner, 
the UT at the conclusion of each integration period was obtained. Hence 
the sidereal time printed out for each integration period was used as 
a time interval measure. The resulting UT's are therefore estimated 
to be accurate to jjD.l sec. 

4. Errors 

Errors in the corrected count data arise from at least three sources. 

(a) The light from the object may have partially missed the spectrometer 
entrance aperture due to seeing variations and guiding/setting errors. 

(b) The effective transmission and seeing of the atmosphere no doubt 
varied to some degree from integration period to integration period 
giving rise to variations in both the object and sky counts. And finally, 

(c) There are inherent variations in the count rate due to photon 
statistics. 
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An estimate of error (a) can be made by examining the first 10 or 
15 minutes of the observations before any darkening occurred. See, for 
example. Part Two, page 141, JII, \A5600 to 5720 Angstroms, and the 
corresponding data in Part Three. The standard deviation of the magni- 
tudes employed to obtain the zero magnitude count is 0.04 magnitudes. 

Since this error will always be a certain fraction of the total counts, 
on a magnitude scale the error will always be a constant. This error is 
immediately visible from an eclipse light curve, and is less than the 
error from sources (b) and (c) in the important refractive tail portion 
of the light curves. It has, therefore, not been included in the error 
bars on the finished eclipse light curves. 

Error sources (b) and (c) may most easily be treated together. 

When the sky is very stable, the standard deviation of the post disappearance, 
uncorrected object count rate scatter about the mean is a good estimate 
of the significance of the uncertainty in the background sky removal. 

It is this error which is shown as error bars on the eclipse light curves 
between the time the telescope was first opened to 200 inches and the 
last point prior to when the curve first fell below la above zero counts. 

Quantitative detail on this error is presented in Part Three for the 
April 5, 1971 eclipse of Europa. During that particular eclipse the 
uncertainty in the background sky due to sources (b) and (c) was only ~1.5 
times the limiting photon statistical uncertainty. 

5. Summary of Data Reduction Procedures 

Table VI-2 summarizes the regions of application of the various data reduc- 
tion procedures. On the first and third events, basically poor data prevented 
use of the preferred curve fit sky removal method. Additional numerical 
detail for the April eclipse of Europa is presented in Part Three. 
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Table VI-2 

SUriMARY OF DATA REDUCTION METHODS 
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BFl . The sky removal method of equation VI-1 was used, employing the balance factors obtained by 
equation VI-2 in the region indicated above. 

BF2. The sky removal method of equation VI-4 was employed. 

Curve Fit: The curve fit method of sky removal was employed. 



VII. DATA ANALYSIS 


Introduction 

Although this report is primarily intended to document only obser- 
vations and basic data reduction, some cursory analysis has been 
performed in this study. The real significance of the observations 
becomes evident when detailed interpretations are made which yield an 
insight into the composition, physical structure, and geometry of the 
Jovian atmosphere and satellites. Although the analysis of the results 
is preliminary, this section has been included as part of the report 
because it demonstrates the potentially significant interpretations 
that can be drawn from this data and from future observations of this 
type. 

A. The Theory of Eclipse Light Curves 

1. The Basic Relationships 

A basic equation for the angle of refraction, ©^, experienced 
by a light ray which tangentially traverses the atmosphere of Jupiter 
in a grazing incidence path is derived by Baum and Code (1953) and is 
given by 

0^ = (2Trroa)''^ (nQ-l)e“®^''l ' ’^o) . (VIl-1) 

(See Table VII-1 for notation.) 

The assumptions implicit in this result are: 

(a) that Jupiter is spherical 

(b) that it is surrounded by an isothermal, homogeneous 
atmosphere 
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TABLE VII-1 NOTATION 


a. . . . Reciprocal scale height, 1/H. [km"^] 

A. ... A constant used in equation (VII-18) to consolidate other 

constants. 

A(x,y) . Fractional beam strength reduction factor physically produced 
by scattering, absorption and cloud top cutoff. A function 
of the wavelength band under consideration and of (x,y) in 
that these define a unique path through the Jovian atmosphere, 
[dimensionless] 

B. . . . A limb darkening constant (See Equation VII-18). Also the 

distance shown in Figure VII-1. [km] 

C. . . . Same as B, 1st usage. 

Dji^. . . Distance from Jupiter to the satellite under study, [km] 

D^j. . . Distance from the sun to Jupiter, [km] 

E. . . . Height above the reference radius, r^, in Jupiter's 
atmosphere in scale height units. E = (r-| - rQ)/H. 
h. . . . Height above the reference altitude, rg. [km] 

H. . . . Scale height at rQ. [km] 

H(h) . . Scale height as a function of h. [km] 

I{x) . . Solar radiance (actually the radiance integrated over z from 
the upper to lower limb at x) per unit of x and per unit of 
solid angle, in the direction of Jupiter, and within the sub- 
ject wavelength band, [watt/sr.km] 

JI-JIV . The Galilean satellites of Jupiter. 
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K-] . . . A consolidation of other constants (See Equation VII-6). 

L(x) . . The fractional solar illumination (relative to the total 

disk) per unit of x, per unit solid angle in the direction 
of Jupiter, and within the subject wavelength band. See 
Equation VII-5. 

n,n(h),n (h). The index of refraction of elevation h and wavelength 

A 

X. [dimensionless] 

n^ . . . The index of refraction at the reference radius r^- Implic- 
tly this is due to H 2 and He and is not near an absorp- 
tion line, [dimensionless] 

N^(h). . The number density of specie x at altitude h. [#/unit 
vol] 

P(h) . . The pressure as a function of altitude h. 

Px(h). . The partial pressure due to specie x as a function of 
altitude h. 

r^ . . . The reference radius at which n = no- [km] 

r] . . . The distance of closest approach of a ray. See Figure VII-1. 

[km] 

Rj . . . The nominal radius of Jupiter. We have taken ro = Rj- [km] 
. . . The radius of the satellite ("moon") under discussion, [km] 

R^ . . . The radius of the sun. [km] 

U. . . . A limb darkening constant, a function of wavelength. See 

discussion after Equation VII-16. [dimensionless] 

V. . . . Similar to U. 
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X. . . . The length axis centered on the sun and shown in Figure VII-1. 

[km] Also a subscript denoting an atmospheric specie, H 2 , 

NH 3 , CH4, He, . . . 

y. . . . The length axis centered on the Jovian satellite orbit and 

shown in Figure VII-1. The satellite need not be centered on 
y = 0 . [km] 

z. . . . The length axis perpendicular to the plane of Figure VII-1. 

[km] 

0 ^ . . . The angle of refraction experienced by a ray due to passage 
through the Jovian atmosphere; positive counterclockwise in 
Figure VII-1. [radians] 

0 ^ . . . The angle at which a ray leaves the sun; positive counterclock- 
wise. See Figure VII-1. [radians] 

A. . . . Wavelength 

Po . . . The atmospheric density at r^. 

<}i/(J)q . . The fractional illumination per unit of y on the y axis 

and on the strip 2 Rj^;| thick in the z-direction in the subject 
wavelength band, [dimensionless] See Figure VII-1. 

0/0Q . . The fractional (i.e., relative to when no absorption, scattering, 
or refraction has occurred) integrated (over entire disk) satel- 
lite brightness in the subject wavelength band; a function of 
satellite position in its orbit about Jupiter, [dimensionless] 
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(c) that a constant effective gravity applies throughout, 
and 

(d) that the effective depth of the atmosphere is much less 
than a Jovian radius. For these conditions the index of refraction is: 

n = 1 + (no - l)e"®^^l ‘ (VII-2) 

Then the further assumption, one well warranted in retrospect, is 

(e) that 0^ is small. 

First, consider only the "diming" of the illumination at the 
satellite orbit due to the spreading of the light caused by differ- 
ential refraction in the Jovian atmosphere. The illumination at y 
per unit of y (see Figure VII-1) and falling in a strip 2R(vi thick 
in the z direction (z is perpendicular to the paper in Figure VII-1) 



Here l(x) is the illumination per unit solid angle and per unit of x 
(i.e., integrated over a strip of the solar disk from 2 = -V^-x2 
to z = ) including the effects of solar limb darkening. 

The factor (2R|v|)/(D5j + Dj|v|) accounts for the spreading, in the z- 
di recti on, of the illumination from the sun. It is just the angle sub- 
tended by the satellite in the z-di recti on as seen from the sun. 

The effect of refraction enters through ay/ae^. 
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ATMOSPHERE 



Parenthetically, meridian concentration in Jupiter's atmosphere 
is neglected and the projected limb of Jupiter on the satellite is 
assumed to be a straight line. The inaccuracy introduced by these 
assumptions is small. 

Since only the fractional illumination is desired, equation 
VI 1-3 is next normalized by its equivalent value when no refraction 
occurs which yields the following relationship: 



I(x), renormalized so as to give the fractional illumination per 
unit of X, is defined as the new function L(x) thus: 


L(x) = I 



(VII-5) 


The analytical form of L(x), including limb darkening, is derived 
in Section 2 below. L(x) will be incorporated into equation VII-4, 
but it is first necessary to evaluate ay/ 30^. 

From figure VII-1 one sees that 


r-| ~ B - (x + D^j 05 ) 


(VII- 6 ) 


Hence, letting K h (ZTir^a)'^ (no " U, equation VII-1 becomes 
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^ .s \ - V - X - Dei 0c - r.) 

OpVX^O^) “ K-^ 6 SJ S 0 


(VII-7) 


and 


Therefore 


y(0£» x) X + (VII-8) 


80 c 


^SJ ® ^SJ ®r ^ °JM 


(VII-9) 


then substitution of VI 1-5 and VII-9 into VI I -4 gives 

Rc 


A. 


f Lix 

I. ^ 


dx 


-p . ^SJ ^JM 

^ hj 


(VII-10) 


Since D^j = 7.78 x 10® km and, for JIV, Dj^ = 1.88 x 10®. to a 
very good approximation 


A. 

i>o 


Ir. 


-Rr 


L(x) dx 
aOj^ 0^(x,y) 


(VII-11) 


The function 0 ^(x,y) is implicitly defined by equations VII-7 and 
VII-8 and in practice, given (x,y), 0 ^ and 0 ^ are found numerically 

by an iterative process. 

Next, absorption, scattering, and cloud top cutoff must be taken 
into account. if their effects are known, incorporation into equa- 
tion VII-11 is a straightforward process. A set of the variables (x,y) 
defines 0 ^ and 0 ^ through (VII-7) and (VII-8). Knowing 0^, r.| 
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can be calculated. Then knowing r-^ the grazing incidence path 
through the atmosphere is defined. If the structure of the density 
variation in the atmosphere is assumed, apart from an overall con- 
stant, the equivalent width expressing the amount of light removed 
from a grazing incidence beam can be pretabulated as a function of 
E = (ri - r^)/!! and the overall number density of each of the 
atmospheric species.' Introducing this extinction factor, A, (E), into 
equation VII-11 gives 


4>0 y 


/ 


L(x)A(x,y)dx 
1 + aDj|^0^(x,y) 


(VII-12) 


To obtain the fractional illumination of the entire satellite 
one merely integrates (VII-12) over the satellite disk obtaining 


J0 

00 


f J L(x)A(x,y) 


This equation gives the fractional illumination in the sense that it 
is the brightness relative to that when no refraction or absorption 
occurs . ■is calculation of this integral for different times and 

different physical assumptions that yields theoretical eclipse light 
curves . 

The weighting factor in the y integral neglects any phase 


function perculiarities on the satellite. A special phase function 
would simply require a different weighting function. 
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2. Effects of Solar Limb Darkening 

L(x) is merely the integral of the surface 


brightness over the solar disk between the 



If the sun were of uniform brightness L(x) would simply be 


L(x) = 



(VII-15) 


However, limb darkening must be taken into account. Allen (1962) 
gives the fractional brightness of any unit area on the solar disk 
relative to that at the center of the disk. 


I\(e)/r,(0) = 1 - U - V +UCOS0 +VC0s2e (VII-16) 

A A 


Here the prime denotes the solar continuum, U and V are functions 

of wavelength tabulated by Allen and 0 is the angle between the 

solar radius vector and the line of sight. In terms of our variables, 

(x,z) I 

I 2 z 

COS 0 = 1 - - ^ (VII-17) 

\ r2 
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Defining three new constants 


A = 1-U-V, B = U, C = V 


We may write 


1\{Q)/I\i0) = A + B COS0 + C COS^e 
From the definition of L(x) it is evident that 


(VII-18) 

(VII-19) 


L(x) = 




y r^(0)/r^(o) dz 


+ R$ + 


(VII-20) 


j* J r^(e)/r^(o) dz dx 


Rs - 


By substitution of 
that: 

2 AR 

L(x) = 


(VII-17) and (VII-19) into (VII-20) one may show 








(VII-21) 


It is immediately seen that (VII-21) has one correct limit. If there 
is no limb darkening, i .e. , A = 1 , B = C = 0, then (VII-21) becomes 


(VII-15). 


In the calculations conducted in the course of the present study 
we have employed Allen's values of U and V interpolating linearily 
between tabulated values to obtain the values at other wavelengths. 
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B. Analysis of the Eclipse of Europa, March 10, 1971 

The theory outlined in Part A has been employed and analysis begun 
of the Europa eclipse observations obtained on March 10, 1971. These 
data were analyzed because they are the best data in existence for 
study of the refractive tail of a Jovian satellite eclipse light curve. 
The analysis presented here is still to be regarded as preliminary. 
Attention is first directed to the continuum channels. 

1. Analysis of Continuum Channels 

Continuum channels are those wavelength bands in which the 
Jovian atmosphere has no absorption features significant at the spectral 
resolution and photometric precision employed. The only sources of 
extinction are Rayleigh and aerosol scattering. Differential refrac- 
tion and cloud top cutoff, together with the physical properties of 
the atmosphere, and these extinction mechanisms, should enable an under- 
standing of the continuum channels. There are enough continuum channels 
so that several may be used to define the parameters. The remaining 
ones may then be used for consistency checks. 

Use of the eclipse light curve technique requires certain data 
describing the Sun-Jupi ter-satel 1 i te system. Summarized here are both 
the parameter values which have been fixed and those which will be 
varied for this analysis. 

Sun: A solar radius of 6.9598 x 10^ km and the wavelength dependent 

solar limb darkening constants, U and V of Allen (1962) have been em- 

O 

ployed. Jupiter has been assumed to be at 7.7836 x 10 km from the 
solar center. 
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Jupiter: As noted previously, an exponential atmosphere is assumed 

in the present analysis. At a reference level in the atmosphere, 71,400 

km (Newburn, 1969) from the center of Jupiter, we have adopted the par- 

20 

tide number density 1.02 x 10 particles/cc, employed by Price and 
Hall (1971). We have also adopted their H2/He ratio of 5.1. These data 
uniquely determine the index of refraction and Rayleigh scattering 
cross section leaving only the scale height, aerosol scattering cross 
section, cloud top altitude and band absorption opacities to be deter- 
mined. 

Satellite: Europa being under study, we have adopted a Jupiter- 

satellite distance of 6.71 x 10^ km (Newburn, 1969). No satellite 
phase function has been employed. Rather, the satellite has been assumed 
to be a uniformly reflecting disk. This assumption, particularly, needs 
further analysis and may account, in part, for the somewhat smaller 
satellite radii yielded by the eclipse technique than those customarily 
employed. In addition, no account of the approximately 11° phase angle 
at the satellite has been made. 

Figure VI I -2 shows how the satellite illumination distribution 
shifts markedly toward the trailing limb of the disk during the first 
five magnitudes of darkening and hence it is evident that the effect 
of the phase function and phase angle may be important in determining 
the satellite radius. During the refractive tail portion the illumin- 
ation is more nearly uniform again. 
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Darkening in Magnitudes 


Illumination distribution 
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Figure VII-2 Satellite Illumination Distribution 


In any case, the satellite radius is the only parameter to be 
determined from the brighter portions of the light curve, but , the 
value determined in this report, or by others using the satellite 
eclipse light curve method, has an inherent uncertainty if consid- 
eration of the phase function is not included. 
a. Satellite Radius 

We begin the analysis by examining the two continuum 
channels, xxl0320-10680 and 6290-6260 A. The data from these two 
channels are reproduced in Figure VII-3. The data are shown down 
to, but not including, the first data point falling below la above 
zero counts. Only the portion of the eclipse light curve brighter 
than 4 magnitudes of darkening for the shorter wavelength channel 
is shown. The slight deviations from a smooth, monotonic decrease 
in the brightest part are felt to stem from the slight loss of light 
outside the spectrometer entrance aperture. The numerical data for 
the subject eclipse show that the sky was stable within ~1.5 times 
the photon statistical limit (see Part Three of this report). The 
wavelength effect on the refractive tail. is quite evident in the 
figure. 

The satellite diameter is determined by the shape of the curves 
during the first 6 magnitudes of darkening. Figure VI I -4 shows a 
superposition of 5 curves, corresponding to different satellite radii 
and \ = 10,500 A. In the calculation of these 5 curves for Europa 
the scale height was taken to be 10 km and the extinction factor. 
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Darkening in Magnitudes 


0 



Time in Minutes 


Figure VI 1-4 Effect of Radius on Theoretical 

Eclipse Light Curve 
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•^ext* equal to unity. However, the scale height and extinction 
do not effect the eclipse light curves until 6 magnitudes of darkening 
It is concluded that the best eclipse radius for Europa is 1200 - 100 
km. This value is used throughout the subsequent discussion. It is 
less than the value of 1460 km cited in Newburn's (1969) review of 
planetary data. 


b. Scale Height ^ Aevosol Extinction and Cloud Top Cutoff 


at XX10320-10680 A 

A single wavelength, three values of the scale height 
and two values of aerosol extinction are used to predict the eclipse 


light curve tails. The values employed were: x = 10500 A; H = 5, 

10, 20 km; f^^^ = 0.5, 1; and no cloud top cutoff. 

The extinction factor, f^^^* has the following meaning: An 

additional scattering cross section, over and above the Rayleigh 
scattering cross section is assumed. Thus, the total cross section is 


ojQj (A,h) = (^jh) + ^ Op (5000 A,h) 


Here, (x,h) = the total scattering cross section, (cm“^) encoun- 

tered by radiation of wavelength x at height h 
above the reference level in the Jovian atmosphere 
0[^ (x,h) = the Rayleigh scattering cross section (cm"^) at 

altitude h and wavelength x. 

Implicitly, this formulation assumes that the ratio of aerosol 
particle number density and gas particle number density is a constant 
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above the cloud tops. This assumption is not obviously true and may 
merit further attention. 

The resulting curves, using these physical parameters are shown 
in Figure VII-5. It is concluded that a scale height of 10 ^2 
an extinction factor of 1.0 best fit the data in the axI 0320- 10680 A 
channel for the location on the Jovian disk sampled by the subject 
eclipse. 

The discrepancy during the latter portions of the refractive tail 
demands attention next. It is postulated that the downturn of the 
refractive tail is caused by cloud top cutoff. Theoretical curves 

o 

with the following parameter were next generated: x = 10,500 A, 

H = 10 km, = 1.0, and cloud tops at 2.5, 2.75 and 3.0 scale 
heights above the atmosphere reference level. A superposition of 
these curves and the eclipse light curve are shown in Figure VII-6. 

From these data we conclude that the cloud tops extend to 26 t 1 
km above the reference level. 

c. Another Continuum Channel: XX6290-6360 ^ 

Next the observed data in the A 6290-6360 wavelength channel 
is fitted. Now, however, only the extinction is varied. Figure VII-7 
show theoretical curves generated with the following parameters: 

X = 6225 A; H = 10 km; cloud tops of 25 and 27.5 km; and = 0.5, 

1.0, and 1.5. It is seen that these data are best fit with a value 
of extinction factor, = 1.5. It is now evident that the continuum 
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Darkening in Magnitudes 



Time in Minutes 


Figure VII-5 Effect of Scale Height and Aerosol Extinction on Theoretical 

Refractive Tails 




Figure VI 1-6 Effect of Cloud Top Cutoff 
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Darkening in Magnitudes 



Figure VII-7 Determination of Aerosol Extinction Factor at X6250 A 
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channels will enable determination of the wavelength dependence of 
the aerosol extinction cross section. 

In subsections a, b, and c., theoretical light curves have been 
fitted to the data in a step by step manner. A satellite radius, 

Jovian atmosphere scale height, aerosol extinction factors, and cloud 
top elevation have all been determined. Figure VII-8 summarizes the 
results. It shows the closest approach layer thickness, central 
elevation, and central ray number density versus time. 

2. Analysis of the Absorption Channels 
Absorption channels are those wavelength bands in which 
detectable CH4, NH^, ... absorption occurs. At the present time, study 
of these absorption data has only begun. However, a few comments are 
indicated. 

An overview of the measurements is shown in Figure VII-9, a three- 
dimensional perspective plot of all the wavelength channels for the 
April eclipse of JII. Table VII-2 lists the wavelengths, anticipated 
features and whether or not observed, as deduced from the 3-d figure. 
Positive identification of a feature is indicated when one wavelength 
channel continually evidences more darkening than adjacent or nearest 
continuum wavelength channels. Shading between the continuum and an 
actual data point on the "cuts" in Figure VII-9 indicates absorption. 

The perspective shown does not immediately convey the significant 
variation, with wavelength, of the continuum height. It is evident 
that the data will yield significant abundance information when processed. 
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Figure VII-8 Refracting Layer Thickness, Altitude and Gas Particle Number Density 
The curves show (left) the thickness of the envelope of rays illuminating 
Europa at closest approach to Jupiter; (center) the altitude of the envelope 
center above the reference level; and (right) the gas particle number density at 
this altitude. The data pertain to the April 6, 1971 eclipse of Europa. 




TABLE VI 1-2 PRELIMINARY COMPARISON OF ANTICIPATED 


AND OBSERVED ABSORPTION FEATURES 

IS AN ABSORPTION 

BANDPASS (ANGSTROMS) ANTICIPATED FEATURES FEATURE OBSERVED? 


3320 

- 

3480 

Continuum 

No 

3640 

- 

3800 

Continuum 

No 

4120 

- 

4280 

Continuum 

No 

4600 

- 

4760 

Continuum 

No 

4920 

- 

5080 

Continuum + CH^ 

* 

5240 

- 

5400 

Continuum + CH 4 

* 

5500 


5560 

NH 3 

Definitely 

5600 

- 

5720 

Continuum 

No 

6290 


6360 

Continuum 

No 

6400 

- 

6440 

H 2 Dipole 

Probably 

6800 

- 

6900 

CH 4 Weak 

No 

7130 

- 

7400 

CH 4 Strong 

Definitely 

7440 

- 

7550 

CH 4 Predicted 

No 

7920 


8035 

NH 3 Predicted 

Definitely 

8160 

- 

8360 

H 2 O Predicted 

No 

9080 

- 

9150 

H 2 S + CH 4 

Definitely 

9240 

- 

9280 

Continuum 

Statistics very bad 
due to narrow bandpass 

9960 

- 

10320 

CH 4 

Very strong! 

10320 

- 

10680 

Continuum 

Possible feature 

10680 

- 

11040 

Continuum + CH 4 

No 


* 


Continuum changing to fast, with wavelength, to tell using only the 
3-d plot. 
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VIII. CONCLUSIONS 

Five Jovian satellite eclipses have been observed with the 200- 
inch Hale telescope and multichannel spectrometer in 20 colors. The 
data has been reduced to eclipse light curve form and statistically 
significant results obtained for up to 10 magnitudes of darkening. 

Markedly different degrees of development of the refractive tails 
were observed from one eclipse to another. 

Analysis of the eclipse light curves obtained during the April 
6, 1971 eclipse of Europa has been initiated. It was found that an 
eclipse radius for Europa of 1200 t 100 km gave excellent agreement 
with the brighter portions of the eclipse light curves. A scale 
height of 10 ^2 *^m, a wavelength dependent aerosol extinction compon- 
ent and cloud tops at 26 t 1 km above the reference level in the 
atmosphere yielded excellent agreement between the theoretical and 

observed light curves. The reference level is that elevation at which 

20 

the gas density equals 1.02 x 10 /cc. 

Ammonia and methane bands and probably the H2, fi rst-order- forbidden , 
pressure induced dipole line at X6420 - 20 K were observed. Hence the 
data inherently contains information which will yield new NH3 and CH^ 
abundance measurements. By analysis of the H2 data obtained, further 
insight into the H2/He ratio may be gained. 

The absolute time accuracy and the observation of both an eclipse 
ingress and egress of Ganymede make available new data on the figure 
of Jupiter. 
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PART THREE 


NUMERICAL ECLIPSE DATA, EUROPA, APRIL 6, 1971 



Introduction 


MECED'MG pag^ 


r r 7 • 


I'^oT nu.m 


Part Three contains a numerical listing of the data for the 
April 6, 1971 eclipse of Europa. The meaning of the various entries 
is explained below. The data for this eclipse is presented because 
it shows beautifully-developed refractive tails. Similar data for 
the other eclipses are available on a limited basis from the authors. 
Explanation of Headings 
(D* TELESCOPE APERTURE RATIO 

The ratio of collecting area after/before opening the telescope 
dust cover to 200-inches. This factor x (UNCORRECTED 
OBJECT COUNT - UNCORRECTED SKY COUNT x BALANCE FACTOR) yields 
the CORRECTED OBJECT COUNT prior to the first data point obtained 
at 200-inch aperture. 

(2) ZERO MAGNITUDE 

The count number used to determine the zero of the DARKENING IN 
MAGNITUDES scale. It is the average of the CORRECTED OBJECT 
COUNT marked in column 3 with a Z, 

(3) FIRST DATA POINT ON 200-INCH 

Indicates at what point in the subsequent data the telescope 
was opened to 200-inch aperture. 

(4) EPHEMERIS ECLIPSE TIME & TIME ZERO 

The date and Universal Time of the eclipse as given in the 


3 



American Ephemeris and Nautical Almanac (1971). The time is 
the zero for the times in the TIME column below. 

BALANCE FACTOR 

The weighting factor relating the sky contribution in the 
object aperture and sky aperture. See Part One, Section VI. 
DOMAIN DETERMINING BALANCE FACTORS 

The region of data analyzed to obtain the BALANCE FACTOR. 
OBJECT MINUS SKY UNCERTAINTY DUE TO PHOTON STATISTICS 
An estimate of the uncertainty in the data due to the photon 
statistics of the sky contributions, i.e.. 

Here 



Cq = average count in the object aperture just after disap- 
pearance of the satellite. C^ = average count in the sky 
aperture at the same TIME. The same region analyzed to obtain 
the BALANCE FACTOR was analyzed to obtain these averages. 

OBJECT MINUS SKY UNCERTAINTY USING CURVE FIT METHOD 
The rms deviation of the actual UNCORRECTED OBJECT COUNT from 
the curve least squares fitted through them in the region bounded 
by and inclusive of the second and third CURVE FIT POINTS numbers. 



This count number is also expressed as a magnitude. Hence 
any magnitude fainter than this differs from zero counts by 
less than la. The 2a and 3a levels above zero counts 
are 0.75 and 1.19 magntidues brighter than the la value. 

(9) SKY REMOVAL METHOD 

Either, CURVE FIT or DUAL APERTURE. See Part One, Section 
VI. This information applied only to the data obtained after 
opening the telescope to 200-inches. The dual aperture method, 
exclusively, was used prior to opening the telescope to 200- 
inches. 

(10) CURVE FIT POINTS 

When using the CURVE FIT method of sky removal the sky in the 
object aperture was approximated by a curve, least-squares-fitted 
to the UNCORRECTED OBJECT COUNT data bounded by the second and 
third CURVE FIT POINTS (e.g., in the first data sheet, .points 
72 to 122). This curve was extrapolated back in time through 
the first point at 200-inch aperture (e.g., in first data sheet, 
point 65). 

(11) POLYNOMIAL COEFFICIENTS 

The constants in the expression SKY BY CURVE FIT “ C-j + C 2 T + 

2 

C^T where T = TIME. If the quadratic curve was found to be 
concave downward, a physically unrealistic occurrence, a straight 
line was least squares fitted through the data and no C^ 

(SECOND ORDER) shows in the printout. 
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(12) NUMBER PLOT, 0 MAG (Columns 1, 2, 3) 

Column 1 gives the data point number, sequentially. Column 2, 

P (or blank) means the data was plotted on the DARKENING IN 
MAGNITUDES versus TIME plots in Part Two. The only points 
omitted (a blank in column 2) were those corresponding to a 
spatial scan or those obtained well after disappearance. 

Column 3, Z (or blank) means the data in column 9 were employed 
in the averaging to obtain the ZERO MAGNITUDE counts. 

(13) TIME (MIN) (Column 4) 

Time interval in mean solar minutes relative to the Ephemeris 
Eclipse Universal Time of disappearance. 

(14) DARKENING IN MAGNITUDES (Column 5) 

The CORRECTED COUNTS relative to the ZERO MAGNITUDE counts, on 
a magnitude scale. 

(15) UNCORRECTED OBJECT COUNT (Column 6) 

The actual count data observed in the object aperture during 
the 10-sec integration period. 

(16) UNCORRECTED SKY COUNT (Column 7) 

The actual count data observed in the sky aperture during the 
10-sec integration period. 

(17) SKY BY CURVE FIT (Column 8) 

In those rows containing an entry this represents the estimate 
of the sky count in the object aperture obtained by the SKY 
REMOVAL METHOD, CURVE FIT. 
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(18) CORRECTED OBJECT COUNTS (Column 9) 

The final number of counts due to the Galilean Satellite with 
the telescope aperture ratio change and background sky corrected 
out. 

(19) CORRECTION FACTOR (Column 10) 

The correction factor used in determining the BALANCE FACTOR. 

See Part One, Section VI'. 


* Refer to first de^ta sheet for clarification of definitions of keyed 
(i.e., (1), (2), ) entries. 
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TELESCOPE APERTURE RATIO = 35.00 BALANCE FACTOR = 0.65989 SKY REMOVAL METHOD, CURVE FIT 
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